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Abstract. In the recent paper 1151 we have introduced a method of studying 
the multi-dimensional Kingman convolutions and their associated stochastic 
processes by embedding them into some multi-dimensional ordinary convo- 
lutions which allows to study multi-dimensional Bessel processes in terms of 
the cooresponding Brownian motions. Our further aim in this paper is to in- 
troduce k-dimensional Kingman- Levy (KL) processes and prove some of their 
fluctuation properties which are analoguous to that of k-symmetric Levy pro- 
cesses. In particular, the Levy-Ito decomposition and the series representation 
of Rosinski type for k-dimensional KL-processes are obtained. 
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1. Introduction, notations and prelimilaries 

The purpose of this paper is to introduce and study the multivariate KL processes 
defined in terms of multicariate Kingman convolutions. To begin with we review the 
following information of the Kingman convolutions and their Cartesian products. 

LetV :— 7'(R+) denote the set of all probability measures (p.m.'s) on the positive 
half-line R+. Put, for each continuous bounded function f on M+, 



(1) / fix)n*i,sv{dx) 



r(g + i) 
V^r(s + i) 

+ 2uxy + y^)^/^)(l - u^-^^^ n{dx)iy{dy)du, 

where fi and v ^ V and 5 = 2(s + 1) > 1 (cf. Kingman[T and Urbanik[19j. The 
convolution algebra (P, *i^s) is the most important example of Urbanik convolution 
algebras (cf. UrbanikflQ]). In language of the Urbanik convolution algebras, the 
characteristic measure, say cTs , of the Kingman convolution has the Rayleigh density 

(2) da,{y) = ^ii±I)!l^y2.+i ^ 1)^2)^^ 

r(s + 1) 

with the characteristic exponent x = 2 and the kernel 

(3) K,{x)^T{s + \)Js{x)l{ll2xy, 
where J six) denotes the Bessel function of the first kind, 

{-lf{x/2Y+'^^ 



(4) Js(x)~S 



fe=0' 



k\T{v + k + l) 
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It is known (cf. Kingman 7J, Theorem 1), that the kernel itself is an ordinary 
characteristic function (ch.f.) of a symmetric p.m., say Fs, defined on the interval 
[-1,1]. Thus, if 6s denotes a random variable (r.v.) with distribution F^ then for 
each t G M+, 

(5) As{t) = EexpiitOs) ^ J cos {tx)dFs{x). 

Suppose that X is a nonnegative r.v. with distribution fi (z V and X is independent 
of 9s . The radial characteristic function (rad.ch.f.) of /i, denoted by fi{t), is defined 

by 

(6) fl{t) = Eexp{itXes) = As{tx)pi{dx), 

Jo 

for every t G The characteristic measure of the Kingman convolution *i^5, 

denoted by as, has the Maxwell density function 

= '^^^(ll\'^' ^''^'exp{-is + l)x'h (0<x<oo). 

and the rad.ch.f. 

(8) asit) =exp{^t^/'i{s + l)}. 

Let P :— V{M.) denote the class of symmetric p.m.'s on M. Putting, for every 
G€V, 



/•OO 

Fs{G) = / F,sG{dc), 
Jo 



we get a continuous homeomorphism from the Kingman convolution algebra {V, *i,s) 
onto the ordinary convolution algebra (V,*) such that 

(9) Fs{GiHjG2} = (FsG,) * {FsG2) {G,,G2eV) 

(10) Fsas = A^(0,2s + 1) 

which shows that every Kingman convolution can be embedded into the ordinary 
convolution 

Denote by M+'^,A; = 1,2,... the k-dimensional nonnegative cone of R*^ and 
'P{M.+'') the class of ah p.m.'s on R^'^ equipped with the weak convergence. In the 
sequel, we will denote the multidimensional vectors and random vectors (r.vec.'s) 
and their distributions by bold face letters. 

For each point z of any set A let Sz denote the Dirac measure (the unit mass) at 
the point z. In particular, if x = (xi, X2, ■ ■ ■ , Xk) G R''"^, then 

(11) Sx ^ X ^ ■ ■ ■ X 6xk, (k times), 

where the sign " x " denotes the Cartesian product of measures. We put, for 
X = (a;i, • • • ,Xk) and y = {yi,y2, ■■■ ,yk) & 

(12) 6xOs,k6y = {5x^Os6y^} ^ {6x20s6y2} X ■ ■ ■ X {(5^ ^ o ^ (5j^ J , {k times), 

here and somewhere below for the sake of simplicity we denote the Kingman con- 
volution operation *i,5,<5 = 2(s -I- 1) > 1 simply by 0^,5 > ^. Since convex com- 
binations of p.m.'s of the form pTjl are dense in P(R"'"'"') the relation p2|l can be 
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extended to arbitrary p.m.'s Gi and G2 G 'P(M+'=). Namely, we put 

(13) Gi 0.,fc G2 = J J 6^ Os,k <5yGi(dx)G2(dy) 

R+fc 

which means that for each continuous bounded function cf) defined on 

(14) J 0(z)Gi 0.,fc G2(dz) = JJ { J 0(z)(5xO.,fc^y(rfz)}Gi(dx)G2(dy). 

In the sequel, the binary operation Os.fc '^iH called the k-times Cartesian prod- 
uct of Kingman convolutions and the pair (7'(M+'"'), Os,fc) will be called the k- 
dimensional Kingman convolution algebra. It is easy to show that the binary op- 
eration Os,fc is continuous in the weak topology which together with ([T]) and (fT3|) 
implies the following theorem. 

Theorem 1. The pair {V{M.'^''),Qs,k) is a commutative topological semigroup with 
So as the unit element. Moreover, the operation Os.fc is distributive w.r.t. convex 
combinations of p.m. 's in 7'(R+'^). 

For every G € 'P(R+'') the k-dimensional rad.ch.f. G(t),t = {ti,t2,- ■ -tk) £ 
is defined by 

(15) G(t)= j nA.(t,x,)G(dx), 

where x = (a;i, X2, • • • Xk) € R"*"*^. Let ©g = 0^,2, • • • , ds,k}^ where 8s,j are inde- 
pendent r.v.'s with the same distribution Fg. Next, assume that X — {Xi, X2, Xk} 
is a k-dimensional r.vec. with distribution G and X is independent of 0^. We put 

(16) [©., X] = 0^,2X2, ■■■,ds,kXk}. 
Then, the following formula is equivalent to (fTS]) (cf. [T3]) 

(17) G(t) = £;e*<*^[®-^l>, (teR+'=^). 

The Reader is referred to Corollary 2.1, Theorems 2.3 & 2.4 [T3] for the principal 
properties of the above rad.ch.f. Given s > — 1/2 define a map Fg^k '■ 7^(R~'''') 
r{R'') by 

(18) i^s.fc(G) = J {T,,Fg) X (T,,F,) X ... X (T,,F,)G(dc), 

WL+k 

here and in the sequel, for a distribution G of a r.vec. X and a real number c we 
denote by TcG the distribution of cX. Let us denote by Vs,kO^'^^) the sub-class of 
7'(R'=) consisted of aU p.m.'s defined by the right-hand side of p^ . By virtue of 
P^ - (fT5|) one can prove the following theorem. 

Theorem 2. The set VsmO^'^'') is closed w.r.t. the weak convergence and the 
ordinary convolution * and the following equation holds 

(19) G{t) = T{F,^k{G)){t) (teR+'=) 

where .?^(K) denotes the ordinary characteristic function (Fourier transform) of a 
p.m. K. Therefore, for any Gi and G2 G R^*' 

(20) Fs,k{Gi) * F,,fc(G2) = Fs,k{Gi Os,k G2) 



4 



THU NGUYEN 



and the map Fs^k commutes with convex combinations of distributions and scale 
changes Tc,c > 0. Moreover, 

(21) F,,fc(I],,fc) = 7V(0,2(s+l)I) 

where Tig.k denotes the k-dimensional Rayleigh distribution and N{0, 2(s + 1)1) is 
the symmetric normal distribution on MJ^ with variance operator R = 2(s + 1)1,1 
being the identity operator. Consequently, Every Kingman convolution algebra 
(V{M.'^''),Q)sk) "is embedded in the ordinary convolution algebra (7^s fc(R+'^), 
and the map Fg^k stands for a homeomorphism. 

Let us denote by £ = {e = (ei, 62, . . . , Cfe), Cj = zL,j = 1,2,..., k}. It is con- 
venient to regard the elements of £ as sign vectors. Denote R^*^ = {[e, x] : x e 
where [e,x] :— (ciXi, 62^2, . . . , e^Xfe). Then the family {R+''',e € £} is a 
partition of the space R*"'. If X is a k-dimensional r.vec. with distribution G, the 
k-symmetrization of G, denoted by G, is defined by 



(22) G^lj^^eG 



2k 

ee£ 

where the operator Se is defined by 

(23) S'e(x) = [e,x] xeR'^ 
and the symbol S^G denotes the image of G under S^. 

Definition 1. We say that a distribution G G 7'(R'^) is k-symmetric, if the equation 
G = G holds. 

Definition 2. A p.m. F € 7^(R+''') is called Q)s,k — infinitely divisible (Q)s,k~ID), 
if for every m—1, 2, . . . there exists F,„ G P(M"'"'^) such that 

(24) F = F,,„ Os,k F,„ Os^k ■ ■ ■ Os,k F„, (to times). 

Definition 3. F is called stable, if for any positive numbers a and b there exists a 
positive number c such that 

(25) T,F Os.k TbF = T,F 

By virtue of Theorem [2] it follows that the following theorem holds. 

Theorem 3. A p.m. G is Qs.fc ~ID, resp. stable if and only if Hs.k{G) is ID, 
resp. stable, in the usual sense. 

The following theorem gives a representation of rad.ch.f.'s of Os.fc— ID distribu- 
tions. The proof is a verbatim reprint of that for ([14], Theorem 2.6). 

Tiieorem 4. A p.m. fi G ^£'(Os.fe) */ '^nd only if there exist a a-finite measure 
M (a Levy's measure) on R"*"*^ with the property that M(0) = 0, M is finite outside 
every neighborhood of and 

(26) / Jl^M(dx) < 00 

7R+fc 1 + ||x||^ 

and for each t — (ti, tk) G R+'^ 

(27) - logA(t) = j^Jl - n A.(t,a;,)}i±M^M(dx), 
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where, at the origin 0, the integrand on the right-hand side of \27y is assumed to 
be 



(28) Zim||,||_.o{l - n ^s{t,x,)}^^ - ^U>^^t 



3 = 1 " 



2 

2 

3 



for nonnegative Xj,j ~ l,2,...,fc. In 'particular, if M — 0, then [i becomes a 
Rayleighian distribution with the rad.ch.f (see definition^ 

k 

(29) -logA(t) = -5]AX', teK+^ 

for some nonnegative Xj,j — 1, fc. Moreover, the representation |i?7| ) is unique. 

An immediate consequence of the above theorem is the foUowing: 

Corollary 1. Each distribution fi € ID(Qs,k) is uniquely determined by the 
pair [M,A], where M is a Levy's measure on R^'' such that M(0) = 0, M is 
finite outsite every neighbourhood of and the condition 1126]) is satisfied and 
A := {Ai, A2, • • • Afc} G M^'^ is a vector of nonnegative numbers appearing in \29(l . 
Consequently, one can write /i = [M,A]. 

In particular, i/M is zero measure then ^ = [A] becomes a Rayleighian p.m. on 
M.^'' as defined as follows: 

Definition 4. A k-dimensional distribution, say f^s,k, is called a Rayleigh distri- 
bution, if 

(30) Ss.fc = Us X Us X • • • X tJs {k times). 

Further, a distribution F G V{M.'^'') is called a Rayleighian distribution if there exist 
nonnegative numbers Xr,r = 1,2 ■ ■ ■ k such that 

(31) F = {Ta.ctJ X {Ta.ctJ X ... X {Tx.as}. 

It is evident that every Rayleigh distribution is a Rayleighian distribution. More- 
over, every Rayleighian distribution is Os,fc— ID. By virtue of ([7]) and (|30p it follows 
that the k-dimensional Rayleigh density is given by 

')k(„ I 1 ]k{s+l) 

(32) g(x) = n)-^i r'^(s + i) ^r'^^p{-(^ + i)IWI'}' 

where x = (a;i, X2, . ■ . , Xk) G R^'^' and the corresponding rad.ch.f. is given by 

(33) S,,fc(t) =i;xp(-|t|V4(s + l)), tGM+'=. 
Finally, the rad.ch.f. of a Rayleighian distribution F on M+'^ is given by 

(34) m = Exp{-\j2>^V3) 

3 = 1 

where Aj, j = 1, 2, . . . , A: are some nonnegative numbers. 
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2. Multivariate Bessel processes 

3. MULTIVRARIATE KINGMAN-LEVY PROCESSES AND THEIR LEVY-ITO 

DECOMPOSITION 

Suppose that fit,t > is continuous semigroup in /Z3(0s,fc), that is for any 
t,s>0 

(35) /it Os,k iJ^s = tJ^t+s 
and {/it} is continuous at i.e. 

limt^ofit = 5q. 

By virtue of Theorem [2] it foUows that {J-'s^kil^-t)} is an ordinary continuous con- 
volution semigroup on R*^. Putting, for each x e and for every Borel subset 
£ of M'=+, 

(36) P(<,£,x) =/itO.,fc<5x(f) 

and using the rad.ch.f. it foUows that the family {P{t,£,x.),t > 0} satisfies the 
Chapman-Kolmogorov equation and, consequently, the formula ([36|) defines tran- 
sition probabilities of a M'^"'"— valued homogeneous strong Markov Feller process 
{Xf,t > 0}, say, such that it is stochastically continuous and has a cadlag version 
(compare [11], Theorem 2.6). 

Definition 5. A M'^'^ -valued stochastic process {Xf,i > 0} is called a Kingman- 
Levy process, if Xf = 

(i) Xo = (P.l); 

(ii) There exists a MJ'~^ — valued homogeneous strong Markov Feller process having 
a cadlag version {Xf,i > 0} with transition probabilities defined by \36\) and Xt = 
X°,i>0; 

4. Fluctuations of Multidimensional Bessel Processes 

Definition 6. Let (Wt,t > 0) be a d-dimensional Brownian motion (d~l, 2, . . .). 
The Euclidean norm of {Wt), denoted by Bt,t > is called a Bessel process. 

It has been proved that Bessel processes inherit the well-known characteristics of 
Brownian motions: They are independent stationary "increments" processes with 
continuous sample paths. The term 'increment' is defined as follows: 

Definition 7. For any s > u the random variable \Ws — Wu\ is called an increments 
of the Bessel process. 

The following theorem gives a Levy-Khinczyn representation of the Bessel pro- 
cess in the sense of the Kingman convolution. 

Theorem 5. The radial characteristic function 4>{x) of the Bessel process [Bt) is 
of the form 

tx^ 

(37) cl>{x) = exp{- } x,t>Q 

4(s -I- 1) 

where d=2(s+l). 
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Since for any s > u the 'increment' of the Bessel process {Bt) is infinitely divisible 
in the ordinary convolution * we have the following representation of the Fourier 
transform of B^-u- 

(38) ^B,,^Ax) = exp{-{s - u)^{x)) 
where tp{^) is a symmetric characteristic exponent 

(39) i}{x) = -(7^ + / (1 - cos xv)U{dv) 



where the measure 11 satisfies the condition 
begin equation 

POO 

(40) / (Tnm(l,a;^)n(da;) < oo. 

Jo 

which implies the following Levy-Ito decomposition. 

Theorem 6. (Levy-Ito decomposition) There exists a Brownian motion X^^^ and 
a compound Poison process X^"^^ independent of X^^^ such that 

(41) Bt = 11^^,11 {t>0). 

Before stating the Wienner-Hopf factorization (WHf) theorem for Bessel pro- 
cesses wc introduce some concepts and notations. The importance of WHf is that 
it gives us information of the ascending and descending ladder processes. We be- 
gin by recalling that for a,(3 > the Laplace exponents K{a, (3) and k{a, 0) of 
the ascending ladder process (£~^, H) and the descending ladder process {L~^,H). 
Further, we define 

Gt = sup{s <t : Xg = Xg} and Gt = sup{s <t : X^ = Xg. 

Theorem 7. (Wienner-Hopf Factorization) Let {Bt,t > 0) be a Bessel process. 
Denote by Bp an independent and exponentially distributed random variable. 

The pairs {Gep,Xep) and (ep — Ge^,Xe^ — Xe^) are independent and infinitely 
divisible, yielding the factorization 

(42) ^^TTTm = ^) ^) ^' ^ ^ ^' 

p — tv + ip{0) 

ip~^,'tjj~ being Fourier transforms and called the Wienner-Hopf factors. 

5. Levy-Ito decomposition of Kingman-Levy processes 
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